Speeding up critical system dynamics through optimized evolution. 
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The number of defects which are generated on crossing a quantum phase transition can be min- 
imized by choosing properly designed time-dependent pulses. In this work we determine what are 
the ultimate limits of this optimization. We discuss under which conditions the production of de- 
fects across the phase transition is vanishing small. Furthermore we show that the minimum time 
required to enter this regime is T ~ 7r/A, where A is the minimum spectral gap, unveiling an inti- 
mate connection between an optimized unitary dynamics and the intrinsic measure of the Hilbert 
space for pure states. Surprisingly, the dynamics is non-adiabatic, this result can be understood by 
assuming a simple two-level dynamics for the many-body system. Finally we classify the possible 
dynamical regimes in terms of the action s — TA. 
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The rapid progress in the experimental realization and 
manipulation of quantum systems [l[ is opening the rich 
and intriguing perspective of the exploitation of quan- 
tum physics to realize quantum technologies like quan- 
tum simulators 0] and quantum computers 0, H| . These 
achievements pave the way to the simulation of con- 
densed matter systems giving the possibility of studying 
different states of matter in controlled experiments 0]. 
Despite the impressive results obtained so far, this is a 
formidable technological and theoretical challenge due to 
the complexity of the systems in analysis and the experi- 
mental requirements. Indeed, the level of control needed 
on the quantum system is unprecedented: one should be 
able to prepare a system in a desired initial state, per- 
form the desired evolution and finally measure the state 
in a very precise way. Moreover, the whole experiment 
should be performed faster than the system decoherence 
time that eventually will destroy any quantum informa- 
tion capability. 

Quantum optimal control (OC) theory, the study of op- 
timization strategies to improve the outcome of a quan- 
tum process, can be an extremely powerful tool to cope 
with these issues It allows not only to opti- 

mize the desired experiment outcome but also to speed 
up the process itself. Traditionally employed in atomic 
and molecular physics [ill, EH , OC has been recently ap- 
plied with success to the optimization of the dynamics 
of many-body systems fl3l - fl5l ]. allowing to achieve the 
ultimate bound imposed by quantum mechanics, the so 
called quantum speed limit (QSL) (To| . Indeed as intu- 
itively suggested by the time-energy uncertainty princi- 
ple, the time required by a state to reach another dis- 
tinguishable state has to be longer than the inverse of 
its energy fluctuations [l7j . This implies that a quan- 
tum system cannot evolve at an arbitrary speed in its 
Hilbert space, but a minimum time is required to per- 
form a transformation between orthogonal states [18l - l22| . 
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FIG. 1: (Color online) Instantaneous excitation energy in the 
LMG model for an optimized (green dashed line, total time 
T ~ Tqsl), a non optimized (red dot-dashed line, T ~ Tqsl) 
and a linear adiabatic process (orange continuous line, T 3> 
Tqsl). Continuous (blue) lines represents the lowest energy 
levels as a function of the driving field F = —t/T. 



For time-independent Hamiltonians this bound has been 
exactly determined [l||; the QSL has been formally gen- 
eralized also to time-dependent Hamiltonians. but so far 
has been computed only in a few simple cases [lj, |23| - 126| . 
A still unexplored, although relevant question is how the 
dynamical crossing of a quantum phase transition (QPT) 
affects this fundamental bound. QPTs indeed represent 
a dramatic change in the low energy sector of a quan- 
tum system and their presence strongly influence its dy- 
namics. The study of the dynamical crossing of phase 
transitions was initially considered in cosmology to inves- 
tigate the formation of the universe [13] ■ More recently 
with the developement of the quantum annealing [28j and 
adiabatic quantum computation]^], a renewed interest 
has been devoted to the subject in condensed matter 
and quantum information [29l - l32| . Here we investigate 
for the first time the QSL of the dynamics of a first order 
QPT in the adiabatic version of Grover's search algo- 
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rithm (GSA) [H H| and of a second order QPT jl| 
in the Lipkin-Meshkov-Glick (LMG) model. Specifically 
we consider the problem of converting the ground state 
on one side of the critical point into the ground state on 
the opposite side in the fastest and most accurate way 
by selecting an optimal time-dependence of the control 
field. We emphasize here that the evolution induced by 
the optimized field is non-adiabatic, as shown in Fig. [TJ 
where the scenario is reproduced for the LMG model, and 
an adiabatic and an optimal evolution are compared: an 
adiabatic strategy (orange continuous line) turns out to 
be effective only for a very large total evolution time T 
(namely T 3> HA , with A being the minimum spectral 
gap, as required by the adiabatic theorem[36j]). When 
the total time is reduced, as realistically required in ex- 
periments in order to preserve the phase coherence, adi- 
abaticity fails, leading to an excited state far from the 
target (red dot-dashed line). However relaxing the con- 
straint of adiabaticity and allowing a more general non 
adiabatic evolution, with OC it is possible to reach the 
desired goal with a fast dynamics (green dashed line). 
Quite surprisingly our study shows that the outcome of 
the dynamical process optimization for the many-body 
systems analyzed is independent from the specific model 
and analogous to that of a two-level system, as sketched 
through the good rescaling of the data in Fig. [5] We in- 
terpret this result as the natural manifestation of the in- 
trinsic metric of the Hilbert space for pure states [23|, [I?} , 
as discussed in Sec. Ill Al Furthermore, studying the QSL 
as a function of the system size, we show that the speed 
up obtained by the adiabatic GSA [H, HU can be re- 
produced and extended to other models with optimized, 
non-adiabatic evolutions. Finally, we introduce the ac- 
tion s = XA as a parameter to characterize the evolution 
of a quantum system and we find that the QSL identifies 
a new dynamical regime, as discussed in Sec. IIIBI and 
summarized in Fig. [5j 



I. MODELS AND OPTIMIZATION 

We study two paradigmatic critical systems, the adia- 
batic GSA [33| and the LMG model [39[ and we compare 
them with the Landau-Zener (LZ) model to better under- 
stand the physics of the process. The GSA Hamiltonian 
is given by: 

H GSA = (l_r(i))(J-|^XW) 

+ r(t)(/-|ite><|iM) (i) 

where the initial state is an equal superposition of all A^ 
basis states \i), i.e. — (J^- \i))/y/N and the final 
target is the specific marked state we want to exctract 
from the database (in our simulations \^g) = 1 10...0) 
without loss of generality). The system undergo a 1st 
order QPT at a critical value of the transverse field T c = 
0.5 (from now on we set h = 1). The gap between the 
ground state and first excited state closes polynomially 
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FIG. 2: (Color online) Infidelity / as a function of the adi- 
mensional scaling variable T/T* for the LMG (red squares), 
Grover (blue circles) and the LZ model (green triangles). 
Data correspond to half of the maximum size analyzed (N = 
64). 

with the size at the critical point: Aqsa ~ A -1 / 2 . 
The LMG Hamiltonian instead is written as: 

N N 

h lmg = _y^ ,/, ; a>; - r(t) £ at, (2) 

i<j i 

where N is the number of spins, erf's (a = x, y, z ) are 
the Pauli matrices on the ith site and Jij — 1 /N (infinite 
range interaction). The system undergo a 2nd order QPT 
from a quantum paramagnct to a quantum fcrromagnct 
at a critical value of the transverse field |r c | = 1. The gap 
between the ground state and first excited state closes 
polynomially with the size at the critical point: Almg ~ 
A~V 3 . We chose as initial state the ground state (GS) at 
I\- 3> 1, i.e. the state in which all the spins are polarized 
along the positive z-axis (paramagnetic phase) . As target 
state has been chosen the GS at T = 0. 
Finally the LZ Hamiltonian that we use as a reference 
model is 

H LZ = T(t)a z +u>a x , (3) 

where the off-diagonal terms give the amplitude of the 
minimum gap Alz = 2w at the anticrossing point T = 0, 
here assumed to be at t = 0, |4(| . In this case the 
initial state is the GS for T(—T/2) = —To and the tar- 
get is the GS for T(T/2) = T , that is —in this effective 
model — to transform the initial GS into the initial ex- 
cited state in an optimal and fastest way. The systems 
analysed have been summarized in the left side of TablcQ] 
For all the models considered our goal is to find the op- 
timal driving control field T(t) to transform the initial 
in the goal state in a given total time T. At the limit 
when the gap closes (the thermodynamical limit for GSA 
and LMG) adiabatic dynamics is forbidden in finite time 
due to the adiabatic condition T ^> A -1 36]: however, 
for finite size systems, an adiabatic strategy might be 
successful. Here we relax the adiabaticity condition, ex- 
ploring a different regime of fast non adiabatic trans- 
formations. Given the total evolution time T, we use 
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TABLE I: Hamiltonian models and predicted scalings in the LZ approximation for the linear and optimal quenches. They 
should be compared with the results of Fig. [4] The functions diverge for the infidelity J — > 0. 




FIG. 3: (Color online) Schematic representation of the ray 
\ifje) evolving along the geodesic connecting the initial state 
\ipi) and the target \i/jg)- 

optimal control through the Krotov's algorithm to find 
the optimal control field T(t) to minimize the infidelity 
I(T) = 1 - \{ip G \iji{T))\ 2 at the end of the evolution, i.e. 
the discrepancy between the final and the goal state @. 
The determination of r opt (t) can be recast in a minimiza- 
tion problem subject to constraints determined by look- 
ing for the stationary points of a functional £[ip, ip, x, L] 
in which the auxiliary states \x(T)) = \4>g)(iPg\' 1 (>(T)) 
play the role of a continuous set of Lagrange multipliers 
to impose the fulfillment of the Schrddinger equation at 
each time during the dynamics, as described in details 

in mm. 

II. RESULTS 

A. Hilbert space metric and optimization 

Previous studies revealed that only when the to- 
tal evolution time exceeds a certain threshold, by iterat- 
ing the algorithm it is possible to reduce arbitrarily the 
value of the final infidelity /. In order to identify such a 
threshold, we fix a target value of the infidelity I* ~ 10~ 3 
and we determined the minimum total evolution time T* 
for which it is possible to satisfy our goal. In Fig. [2] we 
show the value of the infidelity for the optimized pro- 
cess as a function of the rescaled time T/T* for different 
models. The first observation is the presence of a sharp 
threshold, thus, T* can be considered a reliable estimate 
of the QSL for the process considered. The second strik- 
ing feature is the rescaling of the data onto the function 
I = cos 2 (T/T*). We interpret this general behavior as a 
manifestation of the Fubini-Study metric [HI, [37| in the 
Hilbert space. The presence of such a metric for pure 



< 10" 

















•-•GSA lin 

O-O GSA opt 
B-o LMG opt 
■ ■ LMG lin 

a-a LZ opt 
^LZlin 






8 6 8 




20 


100 



N 

FIG. 4: (Color online) The action s* = T*A as a function 
of the size iV (for the LZ model we define an effective size 
iV = A -1 ), where A is the minimum spectral gap and T* the 
time required to reach an infidelity 7* ~ 1CF 3 for the linear 
(full symbols) and optimized (empty symbols) driving field 
for the LMG (red squares), Grover (blue circles) and the LZ 
model (green triangles). 



states can be demonstrated following two independent 
approaches: Wootters in Ref. (37} obtained it just from 
statistical considerations; Anandan and Aharonov in- 
stead in Ref. [23[ derived it generalizing the concept of the 
geometrical Berry's phase to generic non adiabatic evolu- 
tions. According to the Fubini-Study metric the distance 
separating two arbitrary pure states is given by the angle 
between the corresponding rays, 6q/2 = arccos |(V'i|V , G)|- 
For orthogonal states this distance is maximal and given 
by 0g = 7r. The shortest path connecting the states l^j) 
and I^g) is then represented by a geodesic in the ray 
subspace. We can indicate a ray evolving along such a 
geodesic with \tpg), where < 9 < 6q, so that [ Vo) = H'i) 
and \tpe G ) = \ipc), as sketched in Fig. EH For such a ray 
it turns out |(^e|^G)| 2 = cos 2 ((6» G - 0)/2) or, for the 
infidelity, Ig = sm 2 ((0G — &)/2). Numerically can be ver- 
ified that for states on opposite sides of a QPT 0q — tt; 
substituting this value in the expression of the infidelity, 
we obtain 1q ~ cos 2 (0/2). Making the identification 
9 = ttT/T* this last formula is in perfect agreement with 
the data of Fi g. [2] and with the results for a two-level 
system in Ref. [13| . An optimized evolution then can be 
interpreted as a uniform motion along a geodesic with 
speed ir/T*. 
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B. Dynamical regime classification 

In order to establish a classification of dynamical 
regimes and to understand the speed up that optimized 
evolutions gain with respect to non optimized pulses, we 
introduce the action s = TA, obtained through the prod- 
uct of the total evolution time with the minimum spectral 
gap. Notice this is a quite natural way to characterize a 
dynamical process: s » 1 corresponds to a slow evolu- 
tion, for which in principle an adiabatic dynamics could 
be achieved; s < 1 instead characterizes a fast evolution, 
for which adiabaticity is strictly forbidden. In Fig. @] the 
data for optimal driving fields (empty symbols) are com- 
pared with data obtained with a linear time dependence 
(full symbols) T(t) = t/T for the GSA, the LMG and the 
LZ models. We report the product s* = T*A as a func- 
tion of the size N, being A the minimum spectral gap 
and T* the minimal time required to reach an infidelity 
I* ~ 10 -3 . As can be clearly seen in Fig. 21 a linear time- 
dependent T(t) results in an action s* increasing with the 
system size, implying that T* ~ A~ a with a > 1 (full 
symbols). On the contrary, the action s* remains be- 
low the value tt after the optimization (empty symbols). 
Notice that the optimal action is reduced by one to two 
orders of magnitude in Fig. 21 but in general s* 0pt ~ tt 
and s* Lin — > oo for I* — > 0. A simple interpretation of 
the scalings reported in the picture, can be understood 
as follows. For a linear quench, in a first approximation, 
we can assume that the main contribution to the infi- 
delity comes only from the first excited state (29[, and 
the Landau- Zener formula [4^] can be used to give an 
estimate of the excitation probability, i.e. the infidelity: 
I = cxp(— /3A 2 T), with /3 = const, so that by fixing an 
arbitrary (but small) value for the infidelity I* , we have 
T* ~ A" 2 or s* ~ A" 1 . By inserting the gap dependence 
on the size, the scalings reported in Table U are obtained 
for the models considered in this work: they are in almost 
perfect agreement with the numerical data reported in 
Fig. 21 where s* Lin is increasing with the size. The only 
discrepancy is the scaling of the linear LMG model due 
to the fact that the simple linear LZ approximation fails. 
Indeed, here we have that /3 = /3(A) ~ A -1 / 3 , resulting 
in T* ~ A -3 , that is, s* Lin = A 2 / 3 . For the optimized 
process instead, the optimal value s* 0pt = n corresponds 
to a Rabi oscillation between the initial and the target 
state at a frequency Wr = A/2, as clearly shown in Fig. [5] 
or in other words, according to our geometrical interpre- 
tation, the optimal evolution can be seen as a motion 
along a geodesic connecting the initial and the target 
state at a constant speed proportional to A (25j ; an intu- 
itive explanation is provided in the Appendix. As shown 
in Fig. 21 the speed-up we have obtained in our analysis is 
analogous to the speed-up reached through the Grover's 
quantum adiabatic algorithm [33| , from a quadratic to a 
linear dependence on gap of the evolution time T* . In 
the case of the LMG model the gain is stronger, from 
a cubic to a linear dependence, outlining the fact that 
the limit of the optimization is set by a constant value of 
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FIG. 5: (Color online) The different dynamical regimes as a 
function of the action s. 

the action s* = T*A. As a last remark, from the previ- 
ous discussion it can be argued that optimized evolutions 
achieve a substantial speed up only when the minimum 
gap closes polynomially with the size; in the case of an 
exponential closure, even an optimized process leads to a 
total evolution time exponentially diverging with N (4lT - 

We summarize the possible regimes of a quantum 
evolution in Fig. [51 where the final infidelity and the 
timcscalc as a function of the action s are shown. An 
optimized process is characterized by I = cos 2 (s/2) for 
s < 7T and 7 = for s > tt. For a linear non-optimized 
pulse, s < 7r indicates a fast evolution with high defect 
production (I ~ 0(1)) accurately estimated by Kibble- 
Zurek (KZ) theory [29[ or Fermi golden rule (FGR) ap- 
proximation [30| . On the contrary, for s 3> 7r adiabatic- 
ity can be achieved and the infidelity is asymptotically 
vanishing. The QSL then clearly identifies the threshold 
between fast and adiabatic evolutions. 



III. CONCLUSION 

In conclusion, we estimated the time required by an op- 
timized, non-adiabatic process to drive the ground state 
of a many-body system across a quantum phase tran- 
sition. The behavior of the systems analyzed revealed 
to be surprisingly similar to that one of a simple 2-level 
system; we explained the phenomenon through the con- 
nection with the intrinsic geometry of the Hilbert space 
for pure state, interpreting the optimized process as an 
uniform motion along the geodesic connecting the ini- 
tial and target states, as summarized in Fig. [5J This re- 
sult is of particular relevance because establishes a direct 
link between quantum speed limit and Fisher informa- 
tion [22j in the general setup of an evolution driven by 
a time dependent Hamiltonian. We demonstrated that 
the QSL for the dynamical processes analyzed scales as 
the inverse of the critical gap, significantly improving the 
result obtained with a non optimized evolution. Such a 
speed-up is a general feature and can be interpreted as 
the extension of the Grover's algorithm speed-up to the 
models considered. Finally, introducing the action s, we 
provided a classification of the possible QPT dynamics. 
We mention that the understanding of these optimal pro- 
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FIG. 6: (Color online) Instantaneous eigenvalues (blue dashed 
lines) and uncoupled energies of the spin-up, spin-down con- 
figurations (red continuous lines) of the LZ model for a linear 
quench, T ~ t/T. The time t represents the boundary be- 
tween Region I with high transformation rate (|t| < t) and 
Region II where the system is in the instantaneous ground 
state (\t\ > t). 

cess and of the fundamental timescales might be used to 
develop new and more efficient optimization strategies 
also in quantum state preparation. 
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Appendix 

In this Appendix wc provide an intuitive explanation 
for the scaling obtained in Sec. Ill Bl The numerical re- 
sults suggest that the optimized process is equivalent to 
perform a rotation of the initial state into the target at 
a constant angular speed uj = A/2, where A = f(N) 
is the critical gap of the finite size many-body system 



in analysis. In a first approximation we can recast the 
full many-body problem into a 2-lcvcl effective model, 
described with a LZ like Hamiltonian if LZ [F(i)], where 
the spin-up and spin-down states play respectively the 
role of the full many-body initial and target state and 
the tunability of the diagonal element F(i) mimics the 
possibility of selecting the instantanous rotational axis. 
Considering Fig.©, intuitively one can expect that the 
initial state does not change significantly while its en- 
ergy is much larger than the off-diagonal matrix ele- 
ments, while viceversa, an efficient population transfer 
occurs when T(t) < uj. Therefore the total evolution can 
be approximated in two distinct regimes: a first one for 
|r(i)| > uj in which H LZ ~ T(t)a z and the initial state 
corresponds to the instantaneous GS. The second regime, 
for |r(t)| <C uj where the two levels are highly coupled 
and H LZ ~ ujg x . The time t marking the boundary be- 
tween the two regimes is implicitely determined by the 
condition 

r(t")~ Wl (4) 

as shown in Fig.©. Under this approximation, the opti- 
mization problem is now easily solved. It is indeed known 
that the fastest possible transition between two orthog- 
onal states with a fixed overlap uj is obtained through a 
Rabi oscillation, i.e. by applying H = uia x for a time 
T c = it/2uj [1, [23| . The complete rotation is then pos- 
sible if the condition T c ~ t holds, where t is obtained 
solving Eq.(j4|). A linear time-dependence of the field, 
r(t) = t/T, gives t ~ ujT, that is T ~ uj- 2 ~ A~ 2 , im- 
plying s* Lin ~ A -1 . The optimization of the time depen- 
dence of T(t) corresponds to the extension of the region 
of high transition rate, then increasing t: the best possi- 
ble result is clearly t ~ T, that is, the transformation is 
effective during the whole evolution. Finally setting this 
condition, we obtain T ~ n/2uj ~ ttA^ 1 and s* 0pt ~ it. 
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